Restoration of chiral symmetry in the large-A^c limit 
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We study the large- A^c behavior of the critical temperature Tc for chiral symmetry restoration in 
the framework of the Nambu-Jona-Lasinio (NJL) model and the linear cr-model. While in the NJL 
case Tc scales as and is, as expected, of the same order as Aqcd (just as the deconfinement 

1/2 

phase transition), in the cr-model the scaling behavior reads Tc cc Nc ■ We investigate the origin 
of the different scaling behavior and present two improvements of the (T-model: (i) a simple, phe- 
nomenologically motivated temperature dependence of the parameters and (ii) the coupling to the 
Polyakov loop. Both approaches lead to the scaling Tc oc Nc- 
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Although quantum chromodynamics (QCD) is a the- 
ory for quarks and gluons, in the vacuum they are con- 
fined inside hadrons. It is, however, expected that at suf- 
ficiently high temperature and/or density a phase transi- 
tion to a deconfined gas of interacting quarks and gluons 
takes place [U, Moreover, it is also expected that this 
deconfinement phase transition is related to the so-called 
chiral phase transition: chiral symmetry is broken in the 
vacuum and restored in a hot and/or dense medium, see 
Ref. 3 and the lattice simulations of Refs. 

The precise connection between the deconfinement and 
the chiral phase transition is not yet clear. This is also 
due to the fact that both transitions can only be precisely 
defined in limiting situations which arc not realized in na- 
ture. Namely, in the context of pure Yang-Mills theory 
(QCD with infinitely heavy quarks) the Polyakov loop is 
the order parameter for the deconfinement phase transi- 
tion Q : the expectation value of the Polyakov loop van- 
ishes for low T and /x (confined matter) and approaches 
unity in the deconfined phase. On the other hand, in the 
limit of zero quark masses the QCD Lagrangian is invari- 
ant under chiral transformations. The chiral condensate 
{qq) is the order parameter for the chiral phase transi- 
tion: it is nonzero in the vacuum, decreases for increas- 
ing T and /i, and vanishes in the chirally restored phase. 
Nature is somewhere in between: the quark masses are 
neither infinite nor zero. The chiral condensate and the 
Polyakov loop are therefore only approximate order pa- 
rameters. 

Since QCD cannot be directly solved, various meth- 
ods are used to perform explicit calculations. Besides 
the already mentioned lattice simulations, effective mod- 
els containing quark degrees of freedom only, such as the 
NJL model |7HlO{. and purely hadronic models, such as 
the linear a- model [111, [l3 1 have been used to study the 
thermodynamics of QCD. Both approaches cannot de- 
scribe the deconfinement phase transition: the degrees 
of freedom of NJL models are deconfined quarks at all 
temperatures and densities, while linear cr-models feature 
hadronic degrees of freedom in which quarks are always 
confined. However, in both approaches the critical tem- 
perature for the chiral phase transition is in agreement 



with recent lattice studies. 

In order to amend these problems, generalizations of 
the NJL-model have been developed recently, in which 
the Polyakov loop has been coupled to the quarks [H, 
(see also Ref. [1^ , in which -besides quarks- also mesons 
are present). Although confinement is not realized in a 
strict sense because mesons can still decay into quark- 
antiquark pairs [l^, an effective description of confine- 
ment is achieved through the behavior of the Polyakov 
loop at nonzero T and /x. Thus, both deconfinement and 
chiral phase transition can be studied in a unified frame- 
work. 

Another widely used approach to study QCD both in 
the vacuum and in the medium is the so-called large- iVc 
limit [13, 13 in which the number of colors Nc is not 
fixed to 3, but is sent to infinity. When enlarging the 
gauge symmetry of QCD from SU{?>) to SU{Nc) with 
Nc S> 3 one obtains a theory which still contains mesons 
and baryons, but is -although not solvable- substantially 
simpler [l^. Relevant quantities can be expressed as a 
series in A^~", so that it is possible to separate large- 
Nc dominant and large- A^c suppressed terms. This is the 
only known approach to understand some phenomeno- 
logically well-established properties of QCD such as the 
Zweig rule. In order for the large-A'^c Hmit to be con- 
sistent, the following scaling of the QCD coupling gqcD 
must be implemented: 
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In this way quark-antiquark meson masses scale with Nc 
and the amplitude for a A:-leg quark-antiquark interac- 
tion vertex scales as Nc''^~'^^^'^ and thus goes to zero for 
A"c — >■ oo. In particular, decay amplitudes are suppressed 
as 1/^/Nc and therefore quark-antiquark mesons are sta- 
ble and non-interacting in the large-A^c fimit. Thus, at 
nonzero T a non-interacting gas of mesons is realized for 
Nc > 3. 

It is expected that the deconfinement and the chiral 
critical temperatures are independent of the number of 
colors, see Refs. [20l . [2ll | and refs. therein. In fact, they 
should be proportional to the only existing QCD scale. 
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^QCD oc N^. In this work we aim to study this issue in 
detail: namely, we investigate the large-A^c behavior of 
the chiral phase transition in the previously mentioned 
NJL and linear cr-model. Quite remarkably, we find that 
these models behave very differently: while in the NJL 
model the chiral phase transition Tc scales as N^, as ex- 
pected, in the linear cr-model Tc scales as ^/Wl and is 
thus not consistent with the above expectations and with 
the NJL result. Note that these results, although proven 
within the simplest possible versions for both the NJL 
and the linear cr-model {N^ — 2, no vector mesons, etc.), 
are based only on general large- iVc properties and there- 
fore hold also when considering more complicated and 
realistic generalizations. 

The reason for the (inconsistent) behavior of the lin- 
ear cr-models at large-TVc is investigated and two im- 
provements are discussed: (i) a phenomenologically 
based modification of the model, in which (at least) 
one coupling constant becomes explicitly temperature- 
dependent; (ii) the coupling of the cr-model to the 
Polyakov-loop degree of freedom. In both ways the cor- 
rect limit Tc oc iV° is recovered. Thus, it is still possi- 
ble to use chiral hadronic models for studying the chiral 
phase transition, although the present study shows that 
some modifications are needed in order to be consistent 
with the large- A^c limit. 

The NJL-model is a quark-based chiral model 
which has been widely used to study the chiral phase 
transition in the medium. It is based on a chirally sym- 
metric four-quark point-like interaction. The Lagrangian 
in the case Nf = 2 reads as function of Nc- 

Cnjl{Nc) = i^{ii^d^-m,)i> + ^ [{i^^Pf + (V5*75^)]' , 

(2) 

where tp* = {u, d) is the quark spinor, ruq is the bare 
quark mass and G is the coupling constant with dimen- 
sion energy whose A^c-scaling G 3G/Nc (following 
from the relation G oc gQcn) been made explicit. 
The quark develops a constituent mass to* which is pro- 
portional to the chiral condensate (qq): to* = — 2G (qq) . 
In mean-field approximation the effective mass to* as a 
function of T and Nc reads [10] 

l^!!^ + ^r2/V+M dk fc2 2tanh(^^^) 
m*^ Nc\ 2) J„ 7r2 VFT^ ' 

(3) 

where a cutoff A has been introduced in order to regular- 
ize the loop integral. Note that the number Nc of quarks 
running in the loop cancels with the factor 1/Nc from 
the coupling constant, such that the dominant term in 
Eq. ^ is independent of Nc- In the chiral limit to, — 
the critical temperature for chiral symmetry restoration 
Tc is obtained as the temperature at which the effective 
mass TO*, and therefore also the chiral condensate (gg), 
vanish. To leading order in Nc it reads 

T.m^^{l^i-^ oc^r (4) 



When TOg > 0, a crossover is realized and the correspond- 
ing (pseudo-)critical temperature, defined as the point at 
which the first derivative \dm* /dT] is maximal, is also 
A'c-independent. This result, based on general scaling 
arguments, does not change when including the s-quark 
and (axial-)vector degrees of freedom. Moreover, it is also 
unaffected by the 't Hooft terms describing the J7a(1) 
anomaly which is suppressed in the large-TVc limit. We 
thus conclude that in all versions of the NJL model the 
critical temperature for chiral symmetry restoration (sec- 
ond order or crossover) is independent of the number of 
colors. It has therefore, as expected, the same scaling as 
the deconfinement phase transition. 

The linear cr-model is a purely hadronic theory con- 
structed from the requirements of chiral symmetry and 
its spontaneous breaking [ill [3 H^j; of which the 
pions emerge as Goldstone bosons in the chiral limit. 

In order to study the large- iVc behavior of the chiral 
phase transition we consider, as in the NJL model, the 
case Nf = 2 in the chiral limit, in which the Lagrangian 
as function of Nc reads: 

C„[Nc) = \{^,^f + \^.''^'-\^f. (5) 

where <&* — {a, tt) describes the scalar field a and the 
pseudoscalar pion triplet if and where the standard scal- 
ing behaviors \ (suppression of the interaction) 
and /i^ — >■ /i^ (constancy of meson masses) have been 
implemented. For /i^ > a nonzero chiral condensate 
ipo = (p{T = 0) = fiy/Nc/3X ^ Nc/'if^ emerges (/^ is 
the pion decay constant). The tree- level masses for the 
sigma and the pions are to^ — 3A/^ — /x^ , to^ = 0. 

Many investigations of the linear cr-model at nonzero 
T have been performed in the past, see e.g. Refs. (23l - [25| 
and refs. therein. Using the Cornwall- Jackiw-Tomboulis 
(CJT) formalism [2^ in double-bubble approximation 
and excluding the trivial solution, we obtain the gap 
equation [l^ 

N f 
= ifiTf - ^ + 3 y (G, + GO . (6) 

The tadpole integrals over the full propagators Go- and 
Gtt of cr and tt meson read (in the so-called trivial renor- 
malization scheme where vacuum fiuctuations are ne- 
glected) 

J " Jo 27r2yFTMf V ^ / J ' 

where Mi is the effective T-dependent mass of either cr 
meson or pion. At the critical temperature Tc the con- 
densate ip{T) and the masses of a and pion vanish. The 
gap equation (|6]) leads to the expression 

Tc{Nc) = V2U^'^Nl/' ■ (7) 
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For Nc = 3 the critical temperature is, as known, 
Tc — v^/tt [1^, but for an infinite number of colors 
the phase transition does not take place: for Nc ^ 3 
a gas of non-interacting mesons is realized and the me- 
son loops, which are responsible for chiral restoration, 
are suppressed. Note that the scaling behavior of Eq. 
([7]) is not a prerogative of the simple Lagrangian ([5]), but 
holds also in more general hadronic models as those of 
Ref. [i2. 

However, Eq. ([7]) contradicts Eq. This mismatch, 
already noticed in Ref. [13] , is puzzling because both ap- 
proaches contain the same symmetries. Moreover, the 
linear cr-model can be obtained as the hadronized ver- 
sion of the NJL model. However, the hadronization pro- 
cedure should be performed for each temperature T and, 
as a consequence, the coupling constants in the linear 
cr-model should be functions of T. In particular, the chi- 
ral condensate (p{T) of the cr-model should not be larger 
than in the corresponding NJL model. In the following, 
we discuss two improvements of the linear a-model which 
repair the mismatch in the A'c-scaling of Tc- 

We first present a simple and phenomenologically mo- 
tivated modification of the linear cr-model which leads to 
the correct large-iVc results. This consists of replacing 
the parameter //^ with a T-dependent function. 



fi' ( 1 



^0 



(8) 



where the parameter Tg ~ Agp ^ oc Nc introduces a new 
temperature scale. We have implemented the quadratic 
T dependence suggested in Ref. [11]. Inserting Eq. (|5]) 
into Eq. ^ and following the same steps that led to Eq. 
(I7|), the critical temperature now reads 



Tc{Nc) = To 1 + 



2f^Nc 



-1/2 



(9) 



and is independent of Nc in the limit Nc — oo: 
lmiN^^ooTc{Nc) = Tq. It is also clear how the me- 
son tadpoles affect the chiral phase transition. Without 
their contribution, the transition temperature is simply 
Tc = To, like in the large- iV^ limit, and independent of 
Nc- The meson tadpoles are responsible for the term 
oc 3/{Ncf^) in Eq. Q and thus lead to a reduction of 
Tc, Tc < To, for any finite value of Nc- In the case 
Nc = 3, using the numerical value ~ 92.4 MeV and 
setting the temperature scale Tq = Aqcd — 225 MeV, 
the critical temperature Tc is lowered to Tc ~ 113 MeV. 
Interestingly, in the framework of sigma models with 
(axial-)vector mesons, one has to make the replacement 
/tt Z/tt with Z ~ 1.67 [12]. This leads to a critical 
temperature Tc ~ 157 McV, which is remarkably close to 
lattice results [1, Q . With the help of the described mod- 
ification the linear cr-model respects the large-A'c limit 
and is compatible with the NJL model. 

Note that we could have also introduced a T-dependent 
coupling constant A(T) instead of a T-dependent mass 
parameter /i(T)^. As long as A(T) does not depend on 



Nc, our conclusions remain unchanged. Note also that 
the tadpoles in Eq. ^ are natural candidates to induce 
the quadratic T-dependence of Ai(T)^ in Eq. ©. How- 
ever, they are proportional to 1/Nc, while in our case we 
have to require that the loop contributions reponsible for 
the T-dependence in Eq. ([5]) are independent of Nc- Nat- 
ural candidates are, for instance, quark loops, like in the 
NJL model. We therefore expect that the quark-meson 
loop model of Ref. [IB] shows the correct large- A^c scaling 

of Te. 

Our second suggestion for improving the linear cr- 
model is to incorporate the coupling to the Polyakov loop, 
defined as 



l{x) = A^-^Tr 



/ ri/T 
V exp \ igQCD J Ao{ 



, x)dT 



where the trace runs over all color degrees of freedom, V 
stands for path ordering, and Aq{t,x) is the zero com- 
ponent of the gluon field Q. In pure gauge theory 
the expectation of the Polyakov loop 1{T) = {l{x)) is an 
order parameter for the deconfinement phase transition: 
I — Q in the confined phase and Z = 1 in the deconfined 
phase, see the review in Ref. [2^ and refs. therein. 

Following Ref. [lOl (for a similar approach see also Ref. 
[3l| ) we couple the cr-model to the Polyakov loop 



C„-MNc) - C„[Nc)+^\d^l\''T^^V{V)~^ <^^WT^ 



where Ca{Nc) is taken from Eq. ([5|) and the Polyakov 
loop is coupled to the meson fields. Moreover, a kinetic 
term and a potential V{1) for the Polyakov field I have 
been introduced. (Since we are only interested in the 
large-iVj, behavior, the precise form of V{1) is irrelevant 
in the following. Terms of the kind ~ /T$^ could also be 
included jSli] but would not affect the overall A^c-scaling, 
although they might change the order of the phase tran- 
sition.) Applying the CJT formalism in double-bubble 
approximation the gap equation for the condensate V'(T) 
now reads 

= if{Tf - § (m' - h^\l?T') ^ij{Ga + G.) , 

from which the following expression for the critical tem- 
perature Tc is derived: 



Tc = 



h^\l{Tc. 



6A 



(10) 



Assuming that l{Tc) is a constant independent of Nc, we 
again obtain Tc oc A^^ ™ the limit A'c — oo. Detailed 
numerical results represent a task for the future. They 
depend on the form of the Polyakov-loop potential and 
on other parameters of the model. However, the impor- 
tant point here is that it is natural to recover the desired 
large-A'^c limit when the hadronic model is coupled to 
the Polyakov loop. The reason for this is that the chiral 
phase transition is triggered by the Polyakov loop (32j . 
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In conclusion, we have shown that models of NJL-type 
and linear a-type predict a different behavior of as 
function of Nc- In the quark- based NJL model Tc is in- 
dependent of Nc- this result agrees with the expectation 
that Tc scales as Aqcd oc N^, just as the deconfine- 
mcnt phase transition. On the contrary, in linear cr-type 
models a scaling Tc cx -\//Vc is obtained. The different 
scaling originates from the particular mechanism which 
restores chiral symmetry in the two models. In the NJL 
model, quark loops are responsible for chiral symmetry 
restoration, which survive in the large- iVj, limit, while in 
the linear cr-model, meson loops induce the chiral phase 
transition, which disappear in the large- iVc limit. 

This mismatch is, at first sight, even more striking be- 
cause the linear cr-model can be derived from the NJL 
model through a hadronization procedure. However, 
since one should perform this hadronization at each T, 
the parameters of the a-model are actually functions of 
T. Indeed, when making (at least) one parameter of the 



cr-model T-dependent, the expected large- Nc limit can 
be easily recovered. This result, although interesting, 
is based on an ad hoc modification of the cr-model: for 
this reason we have also studied a different approach, in 
which -inspired by Ref. 'st?]- we have coupled the lin- 
ear cr-model to the Polyakov loop. In this way the chiral 
phase transition is induced by the deconfinement phase 
transition and, as a consequence, Tc scales as N°. 

We therefore conclude that the study of the chiral 
phase transition within purely hadronic models is pos- 
sible, although care is needed in order to be in agree- 
ment with the large- iVc limit. Detailed numerical studies 
with more realistic models including (pseudo)scalar and 
(axial-)vector mesons will be presented in the future. 
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